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^ ' We propose to measure the CP phase 5cp , the magnitude of the neutrino mixing 

C**) ■ matrix element \U e 3 \ and the sign of the atmopheric scale mass-squared difference 

0^ ' Amjj! with a superbeam by the joint analysis of two different long baseline neutrino 

oscillation experiments. One is a long baseline experiment (LBL) at 300 km and 
the other is a very long baseline (VLBL) experiment at 2100 km. We take the 
neutrino source to be the approved high intensity proton synchrotron, HIPA. The 
neutrino beam for the LBL is the 2-degree off-axis superbeam and for the VLBL, 
a narrow band superbeam. Taking into account all possible errors, we evaluate the 
' event rates required and the sensitivities that can be attained for the determination 

of <5cp and the sign of Am§i . We arrive at a representative scenario for a reasonably 
' precise probe of this part of the neutrino physics. 



I. INTRODUCTION 



The Super-Kamiokande experiments [1] in the past several years, joined by SNO [2] more 
recently, have given strong indications of neutrino oscillation that are corroborated and con- 
strained by a variety of other experiments. These experiments started a new era in the study of 
neutrino physics and offered the best indication to date of physics beyond the standard model. 
To further probe neutrino physics, there are a number of ongoing and planned neutrino oscilla- 
tion experiments. These experiments promise to give a full description of the phenomenology of 
neutrino mixing. The most attractive experiments among the new generation of neutrino oscilla- 
tion experiments are the long baseline (LBL) experiments. They are performed in the controlled 
environment of traditional experimental high energy physics and expected to allow precision mea- 
surements of the oscillation parameters, including the leptonic CP phase. Notably, the recently 
approved superbeam facility, the High Intensity Proton Accelerator (HIPA) [3] , which can provide 
intensive high energy neutrino beams from its 50 GeV proton synchrotron, offers the possibility 
of even more desirable LBL experiments. So far, the possibility of two LBL experiments using 
the HIPA superbeam have been discussed. One is HIPA to Kamiokande at a baseline length of 
about 300 km [4] known as J2K, and the other is HIPA to a detector located 2100 km away near 
Beijing [5,6] called H2B. It is well-known that there are parameter ambiguities that are generally 
associated with oscillation measurements at a single baseline [7,8]. Measurements at more than 
one baseline can be beneficial [6,9,10]; our previous studies [11,12] showed that the joint analysis 
of the J2K and H2B experiments can offer extra leverages to resolve some of these ambiguities. 
Our results, however, also showed that CP violation effects cannot be determined at 3<7 level 
even with the joint analysis considered in the study, in which no antineutrino beams were used. 

Since the leptonic CP phase and mass-squared difference sign are pertinent information in 
the physics of neutrino mixings, which seems to be very different from that of the quark sector, 



it is necessary to find out how to pin down these neutrino mixing parameters accurately in new 
experiments. It has been widely recognized that the neutrino factory [13] is an ideal facility 
for the study of neutrino mixings. However, because of the technical and budgetary challenges 
faced with building a workable neutrino factory in the near future, and because of the availability 
of a conventional superbeam from HIPA in about five years, it is obviously advisable that we 
explore the full potential of the HIPA superbeams. In this work, we examine in further detail the 
measurement of the CP phase and the mass-squared difference sign from the joint analysis at 
the two long baselines with specific HIPA superbeams and more suitable detectors. The neutrino 
beams are the 2-degree off-axis superbeams for the LBL at 300 km and a narrow band superbeam 
for VLBL at 2100 km. We evaluate the event rates and investigate their sensitivities to the CP 
phase and the sign of Am§ x . Taking into account all possible experimental errors, we find that 
a fairly precise measurement of the CP phase, the sign of the mass-squared difference and the 
mixing angle 9i 3 is possible but requires: (1) the joint analysis at the two baselines; (2) that 
both a beam and a beam are needed at 300 km; (3) that a beam is needed at 2100 km 
if Amjj! is negative; and (4) a significant increase in the statistics at both 300 km and 2100 km. 
We find that sin 2 (2^13) can be probed to very small values, depending on the value of the CP 
phase. 

In Sec. II we describe how the simulations are performed. In Sec. Ill our results are presented. 
A brief discussion and conclusion can be found in Sec. IV. 



II. DESCRIPTION OF SIMULATIONS 



A. Parametrization and inputs 

Our oscillation analyses will be restricted to 3 flavors of active neutrinos. The parameters of 
the system consists of 2 mass-squared differences (MSD), 3 mixing angles and 1 measurable CP 
phase. The unitary mixing matrix in the vacuum is parameterized as usual 

(C12C13 C13S12 §l 3 \ 

-C23S12 - C12S13S23 C12C23 - S12S13S23 C13S23 , (2.1) 
S12S23 - C12C23S13 -C12S23 - C23S12S13 C13C23 / 

where Sjk = sm(9jk), Cjk = cos(6jk), and sjk = sm(9jk)e lS<JP , 6jk defined for j < k are the 
mixing angles of mass eigenstates Vj and v^, and (5qp is the CP phase angle. The three mass 
eigenvalues are denoted as mi, 777,2, and 7773. The two independent MSD are Am^ = m 2 — m 2 
and Anijj = m 2 — m\ . 

The inputs of the mixing angles and MSD's are obtained from solar, atmospheric and reactor 
experiments: 

sin 2 (2<9 12 ) = 0.8, Am^ = 5 x 10~ 5 cV 2 , (2.2) 
sin 2 (20 23 ) = 1-0, |Am^| = 3 x 10~ 3 cV 2 , (2.3) 
sin 2 (26»i 3 ) < 0.1. (2.4) 

Note that the sign of Amf^ is unknown. The currently favored Large Mixing Angle solar solution 
requires Sm^ > 0. 

In LBL experiments the neutrino beam interacts with electrons contained in Earth matter [14] . 
In the present study we use the Preliminary Reference Earth Model [15], generally known as 
PREM, for Earth density profiles [16] and numerically integrate the Schrodinger equation that 
descibes the propagation of the neutrino in matter for the treatment of distance dependent matter 
density. However, we note that there exist more sophisticated approaches to Earth matter effect, 
including both an updated average density profile known as the AK135 [17] and treatments of 
uncertainties of the density profile [8,18,19] which can affected the determination of the CP phase 
angle. 



2 



B. Beams and detectors 

The HIPA 50 GeV proton synchrotron beam calls for a power of 0.77 MW in phase I, to be 
upgraded to 4 MW in phase II. The superbeam provided by HIPA can be a wide band beam 
(WBB), a pulsed narrow band beam (NBB), or an off-axis beam (OAB). The WBB contains 
neutrinos with widely distributed energy. In a NBB the neutrino flux is concentrated in a narrow 
range of energies, with maximum energy -E pea k where the intensity is peaked, and the intensity 
decreases rapidly below -E p0 ak- An OAB also peaks at a certain energy, but has a longer high- 
energy tail than that of a NBB. More details of the various beam profiles can be found in Rcf. [20] . 
For the 2100 km baseline we use a NBB [21] with peak energy of 4 GeV. We have also investigated 
the beams of peak energies of 5, 6 and 8 GeV, but found that the 4 GeV beam gives the best 
results. For 300 km we use the 2-degree OAB (2°-OAB) which has a peak energy at 0.8 GeV [22]. 

The detector at 2100 km is assumed to be a water Cerenkov calorimeter with resistive plate 
chambers [23,5] located in Beijing, tentatively called the Beijing Astrophysics and Neutrino 
Detector (BAND). The size of the detector will be 100 kt at the beginning and can be upgraded 
to a much larger one depending on the physics requirements. The detector at 300 km is initially 
the Kamiokande detector of the present size of 22.5 kt and upgraded later to 450 kt.* 



C. Experimental errors 

For the experimental errors we use 3cr throughout this work. All independent errors, statistical 
and systematic, are added in quadrature. For the statistical error we used 3a Poisson errors as 
described in the appendix of Ref. [24], including a background at the 1% level of the rate of the 
survival channel. For the systematic error we assumed that the background is known at the 2% 
level as given in [25]. To estimate the error due to the uncertainty in the measurement of the 
mixing angle #23, we assumed that sin 2 (2#23) is measured via the survival channel at L = 300 km, 
with the event rate given by N(v^ — > v^) ~ 7Vo(l — sin 2 (20 2 3) sin 2 (A)), where A ~ 7r/2 and N is 
the number of events in the absence of oscillations. Then the statistical uncertainty on sin 2 (2023 ) 
is 

<5(sin 2 (2# 23 )) = VN/N m ^J(l-sin 2 (20 23 )sin 2 (A)). (2.5) 

We then find the variation in the rate in question for a 3a deviation in 9 2 3, and added this in 
quadrature to the other 3a errors described above to obtain the total 3a error. 

D. Scenarios 

We consider three scenarios in the present investigation. The scenarios are summarized in 
Table 1. For Scenario I the first stage involves a 5-year experiment with a water Cerenkov 
detector of 22.5 kt at 300 km with the 2°-OAB beam [22] from HIPA at 0.77 MW. This stage 
is contained in the plan of J2K [4]. The second stage of this scenario has HIPA upgraded to 4 
MW, as discussed in Ref. [5], to deliver a NBB beam to the water Cerenkov detector of 100 
kt at L=2100 km to run for 5 years. 

Scenario II has an upgraded 4 MW HIPA and calls for both and beams. The experiment 
at 2100 km is the same as in Scenario I. For 300 km, however, we assume a much larger water 
Cerenkov detector of 450 kt. It will run for 2 years with 2°-OAB beam, and then 6 years 
with the 2°~OAB D u beam. 



*For a more discussion of the 300 km detector we refer to Ref. [4]. 
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Scenario III is similar to Scenario II, but calls for a much larger water Cerenkov detector at 
2100 km, to run either the or P M beam, for example, for 5000 kt-yr. Whether a v or beam 
is delivered to the 2100 km site depends on the sign of Am^, as will be explained below. 



III. RESULTS OF THE SIMULATION 

A. Strategy 

Let us first describe briefly the strategy of our calculation. There are 6 measurable oscillation 
parameters in the 3-flavor neutrino scheme. There are two mass scales; one is the atmospheric 
scale and the other the solar scale. Existing oscillation experiments have determined that the 
two mass scales are widely separated and therefore sensitive to different L/E v regions. The LBL 
and VLBL we are considering are affected by the solar scale only at next-to-leading order, but 
can be strongly affected by the unknown parameters Sep and #13. Hence we will take Am^ 
and #12 to be the values determined in solar experiments, as given in Eq. (2.2). This leaves 4 
parameters to be determined, but lAm 2 ^ and #23 are known already to a fair degree of accuracy 
from atmospheric neutrino experiments. Therefore obtaining #13, 5qp, and the sign of Am^ is 
the main goal of our calculation. 

The current and soon to be online LBL experiments will determine lAmJ-J and sin 2 (2#23) to 
a better accuracy; Sm^i and sin 2 (2#i 2 ) will be measured more precisely by KamLAND [26,27]. 
The first task of the scenarios we propose is to determine lAm^l and sin 2 (2# 2 3) more accurately, 
using the f M survival probability. This will allow us to have as small an uncertainty as possible in 
the determination of the parameters #13 and 5c p- Let us define N a (L) as the number of charge- 
current events involving the a charged lepton at the baseline L. We assume that iV M (300), which 
depends on the survival probability, is used to determine [Am^l and #23 with as small an 
error as possible. Then in the various scenarios A^(300) cannot be used to determine #13 and 
S C p- * 



B. Scenario I (f M beam only) 

In this scenario we assume that only the beam is employed to run at L=300 km and 
2100 km. Since iV M (300) has already been used to determine lAm 2 ^ and sin 2 (2# 2 3), we are left 
with three types of independent measurements for the determination of #13 , (5cp and the sign of 
Am^: AT e (300), iV e (2100) and A^(2100). The measurements of these three types of events form 
a surface in a three-dimensional space when #13 and Sep are varied in their allowed ranges. The 
angle $13 is constrained by the CHOOZ reactor experiment [28] and the CP phase is completely 
unconstrained. Therefore, we take their ranges to be: sin 2 (20i 3 ) = (0,0.1) and 5qp — (0, 2n). 
Such three-dimensional surfaces, which are tube-like, are displayed in Fig. 1. The upper and 
lower surfaces are for negative and positive Am^, respectively. The closed curves around the 
axes of the tubes are traced out by varying Sqp from to 2ir, while the lines running parallel to 
the axes of the tubes are determined by varying sin 2 (2#i 3 ) from 0.01 to 0.1. For fixed sin 2 (2#i 3 ) 
values, we then obtain the ellipses in Fig. 2. 

When 7V e (300), the number of the v e appearance events at 300 km, is measured, it determines 
a closed curve which is obtained from the three-dimensional surface by a cut at a given value 
on the iV e (300) axis. The value of iV e (300) does not determine #13 directly since 5 is unknown; 
for each of the closed curves we obtain a definite relation between 8 and sin 2 (2#i 3 ) when the 



t When new runs at L=300 km with better statistics are made, the improved 7X^(300) will be used to 
update the values of lAm^l and #23- 
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sign of Amjj is given. We show in Fig. 3 two sets of such relations for each Am^ sign. As 
shown, we choose two extreme values of iV e (300), each of which leads to a range of values for 
sin 2 (26 l i3), depending on the sign of Am^. For positive Am 2 ^, the larger ./V e (300) curve limits 
sin 2 (2#i3) to the range (0.06,0.1), while the smaller one corresponds to sin 2 (26>i3) lying in the 
range (0.006,0.01). Similarly, the ranges of the values of sin 2 (2#i 3 ) for negative Am 2 ,! can be 
read off from Fig. 3. 

We plot in Fig. 4 the two-dimensional curves with fixed A^ e (300). Note that the scale of the 
horizontal axis is logarithmic. If the scale was linear, the curves would be ellipses. An open square 
indicates the point on a curve with <5cp=0°, solid square 90°, open circle 180°, and solid circle 
270°. We also show a representative 3cr error bar for each curve. The assignment of statistical 
and systematic errors has been discussed in the preceding section. The total error is dominated 
mostly by the statistical error. One sees that although the sign of Am^ can be determined at 
the 3tr level, there is no sensitivity to the value of the CP phase. In particular, the error in the 
A^(2100) channel is very large in comparison with the range of variation in the number of events 
when <5cp varies; we will encounter similar situation in the next scenario. 

C. Scenario II (y and v beams) 

By including the P M beam aimed at the detector at 300 km, we have two more types of events, 
i.e., A^ g (300) and iVp(300). So, in addition to the three-dimensional surface in the iV e (300)- 
A r e (2100)-A r M (2100) space shown earlier in Fig. 1, we also have surfaces in the spaces iV e (300)- 
ATg(300)-A^(300) and A e (300)-N g (300)-iV e (2100), as shown in Figs. 5 and 8, respectively. With 
several fixed values of sin 2 (26*13) we obtain the curves shown in Figs. 6 and 9. We plot the 
two-dimensional projections of fixed A*^ (300) =1000 and 10000 in Figs. 7 and 10. 

We found in Scenario I that the 2100 km data with just a source can determine the sign of 
Am 2 ,! at 3cr, but cannot measure the CP phase (see Fig. 4). In contrast, as shown in Fig. 7, the 
300 km data using both a v and v source can determine the CP phase in the ranges (ir/2, 3tt/2) 
or (— 7r/2, 7r/2), but cannot distinguish between the two ranges since the measurement is only 
sensitive to sin S. Furthermore, unless Sep is close to n/2 or 37r/2 the sign of Am^ cannot be 
determined once all of the experimental errors, including the error in the determination of 6*23, 
are taken into account, leaving a four- fold ambiguity. The problem lies in the fact that Ap(300) is 
used; as already noted in Scenario I, survival data provide poor resolution to the CP phase, and 
the matter effect is small at the relatively short distance of 300 km. Six more three-dimensional 
plots which will contain either Ap(300), or AT M (2100), or both, can be made, but they are not 
very useful in the present analysis because they involve the survival data. 

In order to obtain good resolution in the sign of the MSD and to distinguish the two ranges 
of the CP phase as discussed in the preceding section, we have to use data of the electron flavor 
only. Hence we need two experiments with different L/E v ratios and one of them should be a 
VLBL for a good sensitivity to the matter effect. This brings us to the combined analysis of 
A e (300), A g (300) and A e (2100), as shown in Fig. 10. The sign of Am§ x can be easily determined 
if sin 2 (26*i 3 ) is n °t too small and the CP phase can be measured with again the ambiguity 
between the two ranges (7r/2,37r/2) and (— 7r/2, ir/2), as in the case of Fig. 7. The problem lies 
in the fact that the resolution in A e (2100) is poor due to the low number of events, while the 
resolution of the 300 km D e is excellent. So we have to increase the statistics at 2100 km. This 
takes us to Scenario III below. 



D. Scenario III {y^ and beams with increased statistics) 

The situation of Scenario II can be improved if the statistics in A e (2100) are significantly 
increased. This can be achieved by using a larger detector and/or running for a longer period of 
time for the A e (2100) measurement. For Scenario III we set the detector size times the running 
time at £=2100 km to be 10 times larger than that of Scenario II, assuming that the number of 
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events can be straightforwardly scaled up with the detector size. The running at 300 km is the 
same as in Scenario II. The resultant two-dimensional plot is shown Fig. 11. In this scenario, the 
sign of Ara^i can be clearly determined at the 3a level, even for iV e (300)=1000 which corresponds 
to a very small sin 2 (2#i3) lying in the range (0.006,0.01), as indicated by the dotted curves in 
Fig. 11. 

If Am 2 ^ is positive, a reasonably accurate determination of 5c p can be made with no 
sgn(Arn3 1 ) or #13 ambiguity, and the beam is not needed at 2100 km even for very small 
sin 2 (2#i 3 ) in the range (0.006,0.01). This is consistent with the results of Ref. [9], where it was 
found that a — ► v e and — ► v e measurement at short distance and a — > v e measure- 
ment at a long distance could resolve parameter ambiguities for sin 2 (26*13) > 0.005. To see the 
sensitivity more clearly for positive Am^j, we replot the results in Figs. 12 and 13 respectively 
for N e (300) =10000 and 1000. We see that for V e (300) = 10000, which corresponds to larger 
sin 2 (26 i i3) (0.06 — 0.1) as shown in Fig. 3, the CP phase can be determined better than 10° at 
3cx for (5cp small or around 180°. The sensitivity deteriorates slowly when 5c p moves away from 
0° or 180°, and the uncertainty becomes of the order of 25° when (5qp is close to 90° or 270°. 

Even for iV e (300)=1000, which corresponds to very small sin (2^13) in the range of (0.006, 0.01), 
the measurement of the CP phase is still reasonably good. It is interesting to note that the 
sensitivity of the CP measurement near 5cp=0° and 180° for 7V e (300) = 1000 is comparable to 
that of the much higher number of events of N e (300) =10000. Hence, in this scenario, either case 
can establish whether or not CP in the lepton sector is violated if 5cp deviates by than 10° from 
the CP conserving points of <5cp=0° or 180°. 

If Ara^ is negative, the — > v e oscillation is the favorable channel to investigate. Hence once 
it is clear that Am^ is negative (see the next section for a detailed discussion), the beam should 
be delivered to 2100 km to run for 5000 kt-yr. The results, which are the counterparts to Fig. 11, 
are shown in Fig. 14. With positron events at 2100 km the CP phase can be well measured. 
To see the sensitivity more clearly, we replot the results in Figs. 15 and 16 for iV e (300) =10000 
and 1000 respectively. The accuracy of the 5c p measurement for Am^ < using the beam 
is about the same as that of the beam for Am^ > 0, although the distinction between Sep 
in the range (— 7r/2,7r/2) and 5cp in the range (7r/2,37r/2) is not as good for a beam with 
Ara^ < 0. 

We have also done the analysis assuming a v or v NBB of peak energy 5, 6, or 8 GeV is 
delivered to the detector at 2100 km. We found that for these cases the ellipses in Figs. 11 to 
15 are much flatter than for the 4 GeV NBB, so that they do not do as well in resolving the 
degeneracy in 5c p- 

IV. CONCLUSION AND DISCUSSION 

We conclude that with a superbeam, such as that delivered by HIPA, the joint analysis at two 
baselines, of which one is an LBL at 300 km and the other a VLBL at 2100 km, can determine 
the Am^ sign and give a reasonably precise measurement of the CP phase and #13. To achieve 
this, both and beams are needed for the LBL experiment. The survival events — > 
and — ► i/p are generally insensitive to the matter and CP effects. 

The initial HIPA beam with power 0.77 MW will run with exposure 22.5x5 kt-yr at a 
detector at i=300 km to obtain both survival events iV M (300) and appearance events V e (300). 
The former is used to improve the determination of the mixing angle #23 and mass-squared 
difference lAm^l, so as to reduce the uncertainty of these crucial input parameters. The latter 
can show the existence of an appearance signal for sin 2 (26>i 3 ) > 0.006 and find a crude relation 
between 5c p and #13 as shown in Fig. 3. 

A detailed determination of the oscillation parameters will require an upgrade of the HIPA 
beam power to 4 MW. Using our studies in this paper as a guide, we suggest as one possibility 
the following experimental steps using the upgraded HIPA beam: 

Stage 1: Deliver a 4 MW v 2°-OAB to a 450 kt detector at a distance of L=300 km for 2 
years. The survival events AT M (300) are used to determine more precisely the parameters 6*23 and 
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|Am3i|. The appearance events iV e (300) are used to refine the relation between 5c p and 6*13, as 
shown in Fig. 3. 

Stage 2: A 4 MW v NBB with peak energy around 4 GeV is delivered to a detector at £=2100 
km, to run for 100x5 kt-yr. The survival and appearance events iV M (2100) and 7V e (2100) are 
used to determine the sign of Am| 1; with the most sensitivity coming from AT e (2100) (see Fig. 4). 

Stage 3: A 4 MW 2°-OAB is delivered to the 300 km baseline detector for 450x6 kt-yr 
and iV^(300) and iV g (300) are obtained. The data can only determine 6c p and 9 up to a 2-fold 
degeneracy because of the poor separation between 6 and n — 6 in the iVp(300) measurement, as 
demonstrated in Fig. 7. 

Stage 4: A 4 MW (9^) NBB with peak energy around 4 GeV is delivered to the 2100 km 
baseline detector for 1000x5 kt-yr if Am§ x > (Am§ x < 0). Then at 3a, the value of 5c p can 
be determined to about 10° for values close to 0° or 180°, or to about 25° for values close to 
90° or 270°. The distinction between 6cp in the ranges (— 7r/2,7r/2) and (7r/2,37r/2) is better 
for Ato^ > than for Am^ < 0. The Am^ > case is shown in Figs. 12 and 13 and that of 
Am§ x < in Figs. 15 and 16. 

It is apparent from our calculation that in order to obtain enough statistics to provide a 
reasonably precise measurement of #13 and 5c p the total detector size and running time have 
to be sufficiently large. We have not attempted a detailed optimization; rather, we offer our 
calculation as an example for illustration. A search is still required to determine the optimal 
conditions for the measurement. Eventually uncertainties in the Earth matter density along a 
given baseline as well as uncertainties in the solar neutrino oscillation parameters On and Am^ 
must also be taken into account. 
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Tabic 1 Different possible scenarios in joint analyses 





L=300 km 


L=2100 km 




beam 


power 


detector size x riming time 


beam 


power 


detector size x rvming time 




(2°-OAB) 


(MW) 


(kt- year) 


(NBB) 


(MW) 


(kt-year) 


Scenario I 




0.77 


22.5 x 5 




4 


100 x 5 


Scenario II 




4 


450 x 2 




4 


100 x 5 






4 


450 x 6 








Scenario III 




4 


450 x 2 


Vp 


4 


5000 for Am§ x > 






4 


450 x 6 


Vp 


4 


5000 for Am§ x < 



9 



N e (300km) 

15000 



2200 



2100 

N (2100km) 



2000 




N e (2100km) 



FIG. 1. Three-dimensional surface in the events space N e (300) — 7V e (2100) — ^(2100) in Scenario 
II with CP phase Sep varying from to 2ir and sin 2 (2#i3) from 0.01 to 0.1. The lower (upper) one 
is for Am| 2 > (< 0). The surface in Scenario I is obatined by scaling iV e (300) axis by a factor 

0.77(MW)x5(ycar)x22.5(kt) l_ 

4(MW)x2(ycar)x45()(kt) — 42' 
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FIG. 2. Same as Fig. 1, but for fixed sin 2 (26»i 3 ) = 0.02, 0.04, 0.06, 0.08 and 0.1 for the elipses from 
left to right. 
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FIG. 4. 7^(2100) versus iV e (2I00) for fixed iV e (300). The open square, open circle, filled square and 
filled circle denote Sep = 0, it, tv/2 and 3n/2, respectively. 3cr error bars at some points are also plotted. 
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FIG. 5. Three-dimensional surface in the events space iV e (300) — iVe(300) — A?p(300) in Scenario II 
with CP phase Sep varying from to 2ty and sin 2 (2#i3) from 0.01 to 0.1. The upper (lower) one is for 
Am| 2 > (< 0). 
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FIG. 6. Same as Fig. 5, but for fixed sin 2 (26»i 3 ) = 0.02, 0.04, 0.06, 0.08, 0.1 for the elipses from right 
to left. 
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FIG. 7. 7^(300) versus iVe(300) for fixed ^(300) in Scenario II. The open square, open circle, filled 
square and filled circle denote Sep = 0, n, n/2 and 3n/2, respectively. 3a error bars at some points are 
also plotted. The dashed lines denote errors caused by 823 uncertainty. Note that only the upper error 
bar in Np changed since we take 6*23 = n/4 and any deviation from this always moves the Np result in 
the same direction; this is an artifact of choosing maximal mixing as our starting point. 
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FIG. 8. Three-dimensional surface in the events space iV e (300) — 7V e (2100) — iVe(300) in Scenario II 
with CP phase Sep varying from to 2ir and sin 2 (26*13) from 0.01 to 0.1. The right (left) one is for 
Am| 2 > (< 0). 
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FIG. 9. Same as Fig. 8, but for fixed sin 2 (20i 3 ) = 0.01, 0.02, 0.03, ■ • •, 0.1 for the elipses from the 
lower-right corner to the upper-left corner. 
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FIG. 10. 7V g (300) versus 7V e (2100) for fixed iV e (300) in Scenario II. The open square, open circle, filled 
square and filled circle denote Sep = 0, n, n/2 and 3n/2, respectively. 3a error bars at some points are 
also plotted. 
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FIG. 11. Same as Fig. 10, but for Scenario III. 
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FIG. 12. Same as Fig. 11, but for 7V e (300) = 10000 and Smlt > 0. The stars denote 10° step of CP 
phase from to ir/2. 
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FIG. 13. Same as Fig. 11, but for iV e (300) = 1000 and 8rn 2 31 > 0. The stars denote 10° step of CP 
phase from to 7r/2. 
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FIG. 14. Same as Fig. 11, but for iV c -(300) versus iVe(2100). 
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FIG. 15. Same as Fig. 14, but for 7V e (300) = 10000 and Sm^ < 0. The stars denote 10° step of CP 
phase from n/2 to ir. 
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FIG. 16. Same as Fig. 14, but for iV e (300) = 1000 and 5m 2 31 < 0. The stars denote 10° step of CP 
phase from n/2 to ir. 
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